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1 Introduction 


Theories of transport by a random velocity field are used in a number of important 
problems in many fields of science and engineering. Examples range from mass and 
heat transport in geophysical flows m, to combustion and chemical engineering 
[M] , to hydrology [T5] and petroleum engineering [20] ■ In each of these examples one 
finds important physical processes which involve the transport of a passive scalar 
quantity in the presence of a complex velocity field that fluctuates randomly on 
length and time scales comparable to those on which the transport process occurs. 
A central goal of transport theories is to understand the effects produced by such 
random velocity fluctuations on the mean (ensemble averaged) transport. 

A characteristic feature of these theories, which are usually based on perturbation 
methods, are infrared (long wavelength or low frequency) divergences in the terms 
of the pertinent perturbation expansions due to long-range (nonintegrable) spatial 
or temporal correlations in velocity fluctuations. See, e.g., [2] jTgl [22] and references 
therein. Typically, the dominant divergences in lowest order are of diffusion 
type, and correspond to anomalous diffusion in the ensemble averaged transport 
equations, for which the effective diffusion coefficient increases with time and is 
divergent as time t — > oo. So, according to these theories the mean concentration, 
u, of a passive scalar field being advected by a random velocity field with strong, 
long-range correlations satisfies, under appropriate conditions, an equation of the 
form HIT. 2? 

Ut = c(t)u xx + F(u, u x , u xx ), c(t) ~ t p as t —> oo, with p > 0. (1) 

It is our intent here and in [Bj to analyze the long-time behavior of solutions of 
equations of the form ©■ The analysis assumes F = F(u) to be superlinear, in 
the sense that F(u) = 0(u a ) as u — > 0, with a > 1. The inclusion of the nonlinear 
term F(u) in (JT]) accounts for situations in which the scalar field is not conservative, 
meaning that its concentration u undergoes changes due to physical, chemical or 


2 


biological processes. 


We are concerned primarily with the interplay between anomalous diffusion 
(measured in terms of the exponent p ) and nonlinearity (measured in terms of 
the exponent a) in determining the scaling behavior of solutions as they decay 
to zero. In the present paper we analyze the situation where the nonlinearity is 
analytic and “supercritical” (or irrelevant), in the sense that F(u) = ^2j> a o-jU J , 
with a > (p + 3)/(p+l). We show that in this case diffusion is the dominant effect 
in the limit t —> oo, and determines the scaling form of solutions with sufficiently 
localized initial data as they decay to zero: 

u(x,t) ~*~ 7 / 2 ^( t(p +i )/2 ) as f —> oo, ( 2 ) 

where 7 = p + 1 and the function on the right-hand side of © is a self-similar 
solution of equation ID with c(t) = t p and F = 0. Thus, a curious phenomenon, 
universality , is characterized: solutions of many different equations, but differing 
only in the nonlinear term F{u ), the higher-order asymptotics of c(t), or both, 
nevertheless share the same asymptotic scaling behavior, given by a self-similar 
solution of the time-dependent diffusion equation ut = t p u xx . Thus, it can be said 
that such equations belong to the same universality class , in that all of the members 
of this class exhibit the same asymptotic behavior, insofar as the scaling behavior 
of their solutions (with initial data in suitable classes) is the same. 

“Subcritical” (or relevant) and “critical” (or marginal) nonlinearities, namely those 
varying as F(u) ~ u a with a < (p + 3)/(p+l) and a = (jp + 3)/(p+ 1), respectively, 
in the limit u —> 0, are analyzed numerically in j2j. This analysis does not assume 
F(u) to be analytic at u = 0. In marked contrast with the supercritical case, in 
the subcritical case the asymptotic scaling behavior of solutions is strongly affected 
by the nonlinear term F(u). In particular, the decay exponent 7 = 2/(a — 1) 
is determined by the leading-order term (u a ) in F(u), and the function on the 
right-hand side of © is now a self-similar solution of the time-dependent reaction- 
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diffusion equation ut = t p u xx — u a . Thus, the phenomenon of universality is again 
visible. Equations differing only in the higher-order asymptotics of F(u), c(i), or 
both, fall in the same universality class. 

The critical case is peculiar. It marks the crossover from a scaling regime controlled 
(mostly) by diffusion (supercritical case) to a scaling regime strongly influenced, 
and in certain aspects determined, by nonlinearity (subcritical case). Thus, in 
the critical case neither diffusion nor nonlinearity prevails, and the scaling regime 
which is observed bears some features of the supercritical one (same 7 and scaling 
function </)), but acquires an extra logarithmic decay factor, the imprint of the 
critical nonlinearity. 

The following heuristic arguments motivate the critical case, a = (p + 3 )/(p + 1). 
Anticipating that solutions u —» 0 as t —> 00 , it is possible to simplify the analysis 
by taking F(u ) = A u a . This can be thought of as being an approximation near 
u = 0. With this choice, and recognizing that it is the large-time regime which we 
seek to understand, we simplify the equation to 

u t = t p u xx + A 

Under the parabolic scaling x 1—> L^ p+ 1 ^ 2 x, t 1—> Lt, u 1—> ZA p+ 1 )/ 2 m, L > 1 , this 
equation becomes 

u t = t p u xx + ^ 1 +( 1 -“)(p+ 1 )/ 2 a-u“. 

Thus, provided a > (p + 3)/(p + l), as we iterate such scaling transformation we end 
up with an equation where the nonlinear term decreases by a factor £ 1 +( 1 -“)( p + 1 )/ 2 
at each rescaling. Consequently, as the number of rescalings n —> 00 (equivalently, 
t —> 00 ), the linear diffusion term dominates the nonlinear term and we may expect 
solutions to decay as t~(P +1 )/ 2 , the rate determined by the linear diffusion. This 
argument fails when a = (p + 3)/(p + 1), which we dub the critical case. 

We now state precisely the main result of this paper. For this purpose we introduce 
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the spaces 


B q = {f: 


f(w ) S C^E) and ||/|| < oo}, q > 1 , 


with norm ll/ll = sup^gR (1 + M 9 ) ( /H + f'( w ) ) 


and 


S(°°) = {it : E x [l,+oo) —> E j u(-,t) € £>,j for all f > 1 and HrtH^ < oo}, 


where ||w||oo = sup t>1 ||it(- 


Consider the following initial value problem (IVP) 

u t = c(t)u xx + XF(u), t> 1, x G 
u(x, 1) = f(x), itK 

and assumptions: 


( 3 ) 


(HI) f G B q for some q > 1; 

( H2 ) c(t) > 0 for t > 1 and c(t) = t p + o(t p ) as t —> oo, with p > 0; 
(i/3) F(it) = J2j> a analytic at u = 0, with a > (p + 3)/(p + 1). 

We shall prove the following. 


Theorem 1.1 Assume (HI) — (H 3). TTien there exists an e > 0 such that, for 
||/|| < e, one can find B C Bsuch that the IVP /ias a unique solution u £ B 
which satisfies, for some constant A, 


lim 

£—►00 


Vt^u(V¥^-,t) - Af*(-) 


= 0 , 


( 4 ) 


with f*(x) = e 


(P+l) _2 


Our proof relies on the Renormalization Group (RG) approach. RG methods were 
originally introduced, and proved to be very useful, in quantum field theory sna 
and statistical mechanics EH EH]. Their application to the asymptotic analysis 
of deterministic differential equations (both ODEs and PDEs) was initiated and 
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developed by Goldenfeld, Oono and collaborators [121 ; T51 119] . See m and [25] 
for detailed accounts. The mathematical aspects of the method were rigorously 
established by Bricmont, Kupiainen and collaborators HE]- See also mm- 

In the RG approach the long-time behavior of solutions to PDEs is related to 
the existence and stability of fixed points of an appropriate RG transformation. 
The definition of an RG transformation involves two basic steps. The first step is 
the integration (solution) of the PDE over a finite time interval; its purpose is to 
eliminate the “small time” information in the problem (coarse graining). The second 
step is rescaling, to change the time scales in proportion to those eliminated (by 
integration), so that a nominally constant time scale is under study. The iterative 
application of the RG transformation progressively evolves the solution in time 
and at the same time renormalizes the terms of the PDE. Once a proper RG 
transformation has been found for a particular problem, these terms are divided 
into two types: neutral and irrelevant, according to whether their magnitude is 
unchanged or decreases with each RG iteration. The irrelevant terms iterate to 
zero and the dynamics at large times is then controlled by the neutral terms. This 
accounts for the observed universal scaling behavior of solutions as they decay 
to zero. Thus, the RG provides a natural framework in which to understand 
universality. 

Our results contribute to a large body of literature devoted to the study of long-time 
asymptotics of nonlinear PDEs. The work of Barenblatt and collaborators [3] has 
had a major impact in this field of study, specially in elucidating the importance, 
as well as intricacies, of self-similarity in intermediate asymptotics. Our analysis 
follows closely the one in [5]. See also umm- In spirit, our contribution relates 
also to the work summarized in mm- 

Our analysis can easily be extended in a number of directions. These include 
equations in more than one space dimension, nonlinearities involving u as well 
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as derivatives of u (as in equation ©), and nonlinearities with time-dependent 
coefficients. The modifications needed in each case are straightforward. For 
instance, if the nonlinearity is of the form d(t)F(u ), with d(t) ~ t r as t —> oo and 
F(u) ~ u a as u —> 0, the elementary scaling argument presented above suggests 
that the critical exponent is now a = (p + 3 + 2 r)/(p + 1). With this proviso 
Theorem o should also hold in this case. 

The rest of this paper is organized as follows. In Section[2]we establish the existence 
and uniqueness of solutions for problem Q over a finite time interval. In Section [3] 
we employ the RG approach to extend this result to an infinite time interval and 
obtain the long-time asymptotics of solutions, thereby proving Theorem ll.il 

2 Local Existence and Uniqueness 

In this section we prove the existence and uniqueness of solutions for problem ([3]) 
over a finite time interval using a fixed-point argument. In the next section we 
employ the RG iterative procedure [9] to extend this local result over an infinite 
time interval. In the process we obtain upper bounds that lead to the limit ©. 

As a preliminary remark, we note that generically the constants obtained in this 
and in the next sections depend on q > 1, the function c{t) and the coefficients a,j 
of F(-). However, for simplicity of notation we omit this dependence. Also, without 
loss of generality we assume A £ [—1,1] so that the estimates obtained will be valid 
uniformly with respect to A. 

We start with the definition of certain spaces and operators that will be used 
throughout this paper. For L > 1, we introduce the space 

B ^ = {u : R x [1, L] —> R | u(-,t) £ B q for all t £ [1, L]}, (5) 
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with norm || u\\l = sup ||u(-,t)||. Next, let 



™~f{y) d y, 



with s(t) = J* c(v) du, and define the operator 


T(u) =Uf + N(u). 


( 7 ) 


We shall prove that the operator T has a unique fixed point. This is equivalent to 
proving the existence and uniqueness of solutions to IVP © • Specifically, we prove 
existence and uniqueness of solutions in the ball Bf defined below, provided ||/|| is 
sufficiently small. The introduction of Bf is a necessary step since we only stipulate 
the behavior of F(u) near u — 0. So we define 



( 8 ) 


Theorem 2.1 Assume (H1)-(H3) and let L > 1. Then there is an £ > 0 such that 
the IVP has a unique solution u in Bf, if ||/|| < e. 


The result follows immediately from Lemmas 12.11 and 12.21 since the former 


establishes that T maps Bf into itself and the latter that T is a contraction. We 
note that Theorem 12.11 also holds under the weaker assumption a > 2, instead of 
a > (p + 3)/(p+l) as stated in (H 3). However, the latter is necessary for the proof 
of Theorem II.II in the next section. 


Lemma 2.1 Let L > 1. There is an e' > 0 such that ||^V(u)||i < ||/|| for all 
uG B f , if f G B q and \\f\\ < s'. 

Lemma 2.2 For each L > 1 there is an e" > 0 such that ||JV(u) — 7 V(u)||l < 
||| u — v\\l for all u, v G Bf, if f G B q and ||/|| < e". 
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Before proving Lemmas 12.11 and 12.21 we make some remarks and obtain estimates 
which are needed in the proof and also in Section [3] 


Let us assume that the Taylor expansion of F(u ) at u = 0 has a finite radius of 
convergence p (the case p = oo is less restrictive). We argue that for N(u) and 
T(u) to be well defined it suffices to require that u £ Bf with f £ B q such that 
||/|| < [2(7 g (l + y/s(L))}~ 1 p. Indeed, we have to check that \u(x, t)\ < p for all x £ K 
and t £ [1 , L] and this will follow by comparison of different norms. First, notice 
that B q C L 1 (R) (~1 L 2 ( R) D by elementary Fourier transform calculations. 

For instance, for all x £ R 

\h(x)\ < sup f \h{w)e' wx \dw < d w = C q \\h\\. 

x ^ J-oo 2n Jr 1 + \w\ q 

So, if u g B^ L \ we have it(-, t) £ B q and 

sup \u(x,t)\ < Cq 11U11 

X 

for all t £ [1, L\. Now for f £ B q and u £ Bf we have 

|HU<2(l + v^(I))||/||. (9) 


Define (p(r) = s(t) — s(t — r) (recall s(t) = f±c(v) dv) and J(w,t ) = 
Z *" 1 wip(T)e~ v ^ w2 dr for t > 1 and isei The following estimate holds: 

\J(w,t)\<(t-l)y/s(tj. (10) 

Indeed, taking x = y/ip(r)w and using that xe~ x2 < 1 for all ieR, 


| = 


s/ \p(r)xe x dr 


< 


vM'O dr. 


Since the integrand is a continuous function in [0, t — 1] and ip(r) is an increasing 
function (notice that c(t) > 0 ), we can bound the last integral by — 1 ) = 


(t — l)y/s(t) — s(l) and, since s(l) = 0 , we have (flfil) . 
Also, if q > 1 and w £ R, 

I(w) = j 

J R 


1 + lirl 9 1 + lx — w\ q 


dx < 


C 


1 + \w\F 


( 11 ) 
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where 


C = C(q) = (2 q+1 + 3) / —dx. (12) 

Jr 1+ \x\ q 

Now, motivated by hypothesis (H 3), let C be given by IT2l) and define the sums 


s °(*> = £ Ur 


j>a 


C\ 


3 -1 


a,- z- 


Si ( z ) = 


3 >« 


27T / 

C\ 


']> 


3-1 


3 >“ 


27t / 

3-1 


3-2 


S '2^) = 5Z(^) J\ a j\ zj 2 


(13) 

(14) 

(15) 


Notice that the radius of convergence of these sums is ( 2np)/C < p. We will now 
consider only those functions u such that |u(x,t)| < tt p/C = po for all x € M and 
f £ [1, L\. Invoking (O, it thus suffices to take / such that 


< 


2 C g (l + \J s(L)) 


Po- 


(16) 


Proof of Lemma l2.lt Taking the Fourier transform of N(u) yields 

_ nt-1 ^ 

N(u)[w,t) = a55°3 / e -v , ( T )" J ui(w, t — t) clr. 

3 >« Jo 

Writing ui as convolutions of u, each term in the sum above is of the form 

7 tTyTT [ dre _v(r ^ 2 f it(w-p 1 )u(p 1 -p 2 ) ■ ■-u{pj-i) dp ± -■-dpj-x, (17) 
(2tt) j Jo J Ri-i 

where we have omitted the dependence of u on t — r. Since the absolute value of it 
is bounded by ||u||l/(1 + |w| 9 ), (fTTl) can be upper-bounded by 

1 1 


rt-1 

\3 I d T e ~‘ pl ' T ' ,w2 


dpi ■ ■ ■ dpj-i. 


(2tt)J o 7 kj-i 1 + \w ~Pi\ q l + \pj-l\ q 

Here the integrals over K no longer depend on r and we can bound the exponential 
by one to obtain t— 1 as an upper bound for the integral with respect to r. Therefore, 
using dill, 

--- i — 1 

(18) 


N(u){w,t) < |A| So(Ml), 


1 + \w\ q 
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where Sq is the sum given by m- Similarly, the derivative of G3 with respect to 
w can be bounded above by 


(\2J(w,t)\ +t- 1) 


1 


(2 7 r)J- 11 


1 + \w-pi\ 


i + Ip. 


' j -n 


dpi••• dpj-i 


and using cm) and (fill) we conclude that 


|A^}W)I < |A| 1} (* ~ l)g 0 (||«IU). (19) 


1 + \w\ q 


Estimates (11811 and ([T9]l . together with the monotonicity of s(t) and the inequality 
a > (p + 3)/(p+l), then imply that 

\mu)\\ L < 2\X\{^) + l){L-l)\\u\\lY J (^] 3 'lailNli" 2 . 

j>a ' ' 

We can now bound the sum above by its value when ||u||_l = Po and use (J9|) to 
obtain 

||tf(u)lk<C"|A| ll/ll 2 , (20) 


where 

C' = C\L, q, F,c) = 8 (v^(I) + l ) 3 (L - l)Si(po). (21) 

Finally, recalling that |A| < 1, invoking (TT51) . and defining 

e’ = min jc" -1 , [2C q (l + \/s(L))] _ Vo} , 

where C is given by equation (EH, we conclude that ||AT(u)||l < ||/|| whenever 

ll/ll <£'• ■ 


Proof of Lemma 12.21 Consider functions u and v such that ||u||l < po and 
||"c ||l < Po- Then, 


[N(u) - N(v)]{w,t) = A^aj 

j>a 



dre [t/J — u J ](w, t — t) = A ^ 

j>a 


where Dj can be written as 


pt-i 

Dj = Qj / dre _v, ^’ r ^“ [(u * ■ • • * u) — (v * ■ ■ ■ * v)](w, r). 

(2ny 1 J 0 
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Here there are j — 1 convolutions of u and j — 1 of v. We add and subtract in the 
integrand the term v * u * • • • * with j — 2 convolutions of to get 

nt-i 

Dj = V 7 / e -v( T ) w [(fl — f)) * u * • • • * u)(w, r) dr + 

i. 27T r Jo 

a f t_1 2 

+ 3 / e ~V(T)w [{)*(£*.... ifU-V*---* t ) dT. 

Jo 

The first integral can be bounded, in a manner similar to what was done in the proof 
of Lemma l2Jl by (t — 1)(1 + |w;| ,? )“ 1 ( 27 r) 1_J ( 7 f- 1 | aj ||| M ||Z- 1 || u _ ^o estimate 
the second integral we rewrite it, after adding and subtracting appropriate terms, 
as a sum of two integrals, one of which can be bounded as above and the other can 
be split into two other integrals. This procedure ends after j — 1 steps, when we 
obtain 


Dj< 


(*-l) 

1 + \w\ q 


c ^ 1 M 

(27r)-i -1 




Note that since the norms of u and v in B ( - L '> are less than po, the sum over j > a of 
the right-hand side of the inequality above is convergent. In addition, we can factor 
||u||l or ||v||l and the remaining sum will still be convergent. Similarly, each term 
of the derivative with respect to w of the difference N(u) — N(v) can be written 
as a sum of two integrals, which we bound using the same procedure as before. 
Therefore, 


\\N(u)-N(v)\\ l < C"\\\\\f\\\\u-v\\ L , 


where 

C" = C"(T,g,F,c) =4( V / i(I)+ 1) 2 (L- l)5 2 (po). ( 22 ) 

Since |A| < 1, defining 

e" = min {(2C")- 1 , [2C q {\ + V^))]“Vo} , 
where C" is given by (| 22 |> . the lemma is proved if we take ||/|| < e". I 


We note for future use that, since 5 2 (po) > Si(po), we may take 

C 0 = 8(^iL) + l) 3 (L-l)S 2 (po), (23) 
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and it is enough to consider e in Theorem 12.II defined by 



(24) 


3 Global Existence, Uniqueness and Asymptotic Behavior 

It follows from Theorem 12.II that, given L > 1, there exists an e > 0 such that the 
IVP 0 has a unique solution u in Bf for any / G B q with ||/|| < e. Therefore, 



(25) 


is a well defined element of B q . The right-hand side of (1^51) defines an operator, 
Rl, o, acting on the ball {/ £ B q : ||/|| < e}, which maps the initial condition / to 
/i. We dub Rl,o the Renormalization Group operator associated to problem ([5]). 

The RG operator just defined was introduced in [5]; its iteration comprises the RG 
method for the asymptotic analysis of solutions. The basic idea of this method is 
to reduce the long-time-asymptotics problem to the analysis of a sequence of finite- 
time problems obtained by iterating the RG operator. In more detail, first consider 
problem (J3J) and, as in Section [21 restrict the initial data so that this problem has 
a unique solution. Then, apply Rl ,o to the initial data / to produce /i, the initial 
data for a new, renormalized IVP. It is expected that this procedure can be iterated 
ad infinitum to generate a sequence of finite-time IVPs, whose initial conditions /„ 
are obtained by iterating the RG operator n times. 

We now outline the RG method for the nonlinear problem 0- See also mm- 
Assume that the solution u to IVP 0 is globally well defined and let L > 1 be 
fixed. We consider a sequence {u n }%L 0 of rescaled functions defined by 


u n (x,t) = L n(p+1)/2 M L n t\ , 


(26) 
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where c n (t) = L~ n P C {L n t), X n = £™[p+3-<*(p+i)]/2 A) 


F n (v) = [L n{a ~ j)(p+1)/2 aj 

j><* 


V 


3 


and 


/„(*) = U n (x, 1) = L n(p+1)/2 M (l^p+V^x, L n ) . (28) 


Comparing (1751) and ((4j) , it becomes clear that proving the asymptotic limit may 
be reduced to proving the convergence of {/„}, which motivates the definition of 
the RG operator. Let g € B q and for a given n > 0 assume that the IVP (1271) with 
initial condition g has a unique solution u n . Then, rescale u n {-, L) to obtain 


(Rl,u9){x) = L (p+1)/ V l) , (29) 


which defines the RG operator. The index n in the above definition is justified since 
the operator depends on the evolution equation considered. Now if we consider 
IVP (1771) with initial data f n given by (1751) , then it is an immediate consequence of 
these definitions that the sequence {/„} satisfies 


fo — ui', 1) and ,/n+l — RL,nfn- 


(30) 


Our goal from now on is to make the above heuristic argument rigorous. We shall 
prove that under hypotheses (H 1) — (H 3), if the initial data is sufficiently small, 
problem (1771) has a unique solution for each n so that the iterative RG method can 
be applied to furnish the asymptotic behavior of the solution to IVP ([3]). 


In Lemma 13.11 we proceed as in the proof of Theorem 12.11 to obtain local existence 
and uniqueness of solutions for each problem (1771) . To state the lemma, consider 
the space B <fj ^ defined by © and, if f n is the initial data of problem m, 
define the space Bf n = {u n € B^ : \\u n — u/„|| < ||/ n ||} and the operator 
T n {u n ) ~ Uf n + N n {u n ), where Uf n is the solution of (1771) with \ n = 0 and 


NniUn) {x , t) — A n 



(x-y ) 2 

e”4l« n (t)—n(t-r)j 

\J 47r[s n (f) - s n (t - r)] 


F n (u n {y,t- t)) dydr, 


(31) 
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where 


(32) 


r f fp +1 _ i 

Sn(t) = / Cn(v) dv = -—-h r n (t). 

J1 p +1 

Define also the constant C n by 

C n = 8 (Vsn(L) + lf(L - l)S 2 (p 0 ), (33) 

where ^(po) is given by (fl5l) . with z = po. 

Lemma 3.1 Given n £ N and L > 1, there exists an £ n > 0 such that if \\f n \\ < s n , 
then the IVP \2 7| ) has a unique solution u n (x,t) in Bf n . Furthermore, f n +i given 
by \30 1) is a well defined element of the B q space. 

Proof: (Notice that if n = 0 and /o = /, Lemma IQ1 is just Theorem l2.in We must 
prove that the operator T n is a contraction in Bf n , therefore obtaining a unique 
solution u n in Bf n . First, following closely the arguments in the proof of Lemma liOl 
the constraint L > 1 and the definitions of F n and s n (t) imply that 

||JVn(«n)||i < C n L n ^+ 3 - a ^ +1 ^ 2 \\f n \\ 2 (34) 

and that 

H-^^n) - N n (v n )\\ < C'„i" b+3 -“ (P+1)I/2 ||/r l ||||Un - W„||, 

where C n is given by (1331) . The condition for u n to be in the region of analyticity of 
F n is now that ||/ n || < [2C g (l + i/s„(L))] _ 1 p 0 - Since p+ 3 — a{p+ 1) < 0, defining 

e n = min | (2C„) _1 , [2C q (l + a/s„(L))] _ 1 po} , (35) 

if ||/„|| < e n , we obtain \\N n (u n )\\ L < \\fn\\ and || N n {u n ) - iV n (i;„)|| < \\\u n -v n \\ L 
for all u n ,v n £ Bf n . This proves that the IVP (1771) has a unique solution u n (x,t) 
in Bf n and, therefore, /„+1 = (L^ p+ 1 i/ 2 x, L) is well defined. I 

We have proved that if ||/„|| < e n , then R,L. n fn is well defined. To simplify the 
notation, let v n {x) = N n (u n )(x,L) (cf. (13T1) . where Nq(u) = N(u)). Then, the 
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solution to the IVP (1271) at time t = L can be written as u n (x, L) = Uf n (x, L)+v n (x) 
and we have 

(■ RL,nfn){x ) = R° L fn(x) + L (p+1)/2 V n (L {p+1) /2 x) , (36) 

where R° L = 0 and (R i f n f n )(x) = L^ p+1 ^ 2 Uf n (z/ p+1 )/ 2 :r, L). We see that (l36l) 

splits the RG operator into two parts, which we dub the linear and the nonlinear 
parts. Our analysis focus first on the linear part; the nonlinear part is driven to 
zero under hypothesis (H 3) and, thereby, does not contribute to the asymptotic 
regime, as we shall prove. 

It follows from the definition of R° L n and from the integral representation of uf n 
that if g £ B q is the initial data of IVP (1271) with X n = 0, then the Fourier Transform 
of R °L,n9 is given by 

m°L, n 9)M =g{L-( p +Ww)e- w2 ^/ LV+1 , (37) 

where s n is defined in (1321) . Applying equation (1371) inductively and using that 
s 0 (L n ) + £ n (p+ 1 ) s n (L) = so(L n+1 ) for all n = 1,2,..., it is easy to prove that the 
linear RG operator has the semi-group property. Also, if g = /*, it follows from 
equation m and definition (1571) with n = 0 that 

R(R° L f;)(w) =e -™ 2 /(P+l) e -™ 2 r(L)/Lr+^ ( 38 ) 

where r{L) = r 0 (L) (see (1571) for the definition of r 0 (L)) and we have used that 
fp(w ) = e~ w2Rp+1 \ Taking the inverse Fourier Transform on both sides of 
equation (1551) we conclude that /*(x) is a fixed point of the RG operator if and 
only if r(L) = 0, which is valid if we take c(t) = t p . The next lemma shows that if 
r(t) 0, /* is still the long-time asymptotic limit of R° L f* and, therefore, it will be 
an asymptotic fixed point of the linear RG operator. 

Lemma 3.2 There is a constant M = M(p, q) and an no £ N such that 

\\R%f; - / P 1 < M\L- n ^r(L n )\ (39) 
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for all n > uq. In particular, \\R° L n.f* — f *\\— >0 as n 


oo. 


Proof: Since r(t) = o(t p+1 ), from the semi-group property of Rff and from 
equation (1551) . with L replaced by L n , we have pointwise convergence in Fourier 
space. To prove the theorem, we have to show convergence in B q . This will be done 
at the same time that we estimate the rate of convergence. From (l38l) . 


r(L n ) 


2/1 (P+ 1 ) 


n(p+1) 


and 


-o'Hl <2 


\m 


p + l 


\W\ 


r(L r ‘ 


L n (P+ 1 ) 


i I K*-") 

P+ 1 L n{p+ 1) 


Since r(t) = o(t p+1 ), there exists an no >0 such that |r(L")L _ra ( p+1 ) | < [2(p+l)] _1 
for all n > riQ. Multiplying the inequalities above by (1 + |u;| 9 ) and defining 
M = max UJ (l-|-|w| 9 )e _ [ 2 ^ p+1 ^ ±w2 [2|w;| +w 2 + 2\w\ z j(p+ 1)] we obtain (|39l) . Letting 
n —> oo finishes the proof. I 


In the next lemma we prove that if L is sufficiently large, then the linear RG 
operator Rff n is a contraction in the space of functions g € B q such that y(0) = 0. 
This result will be used to obtain the estimates of the Renormalization Lemma EfTl 
which will allow us to prove Theorem ll.il 


Lemma 3.3 (Contraction Lemma) There exist constants L i > 1 and C > 0 
such that the inequality 

\K, n 9\\ < L- (p+1 )/2 C|| 5 || (n = 0,1,...) (40) 

holds if L> L\ and g £ B q satisfies g(0) = 0. 


Proof: We bound (l37l) and its derivative by 

s n (L) 


! (L) 

LP+ 1 


7- P±1 

L 2 


+ 2 u> 


LP+ 1 


9 


r £+1 

L 2 


+ L~(p+1)/2 


9 


r £±i 

L 2 


From (13^1) . it is easy to see that r n (t) = L ”( p + 1 ) [r(L"f) — r(L n )]. Now, since 
r(t) = o(t p+l ), we conclude that there exists an f o > 1 such that, if L > Lq, then 
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| L ^ P+1 V„(L)| < [2(p + 1)] 1 for all n. This together with the definition of s n 


furnishes 


< s n {L) < 


(n = 0 , 1 , 2 ,...) 


(41) 


6 (p + 1) ^ LP+ 1 - 2 (p + 1) 
if L > Lj = max{i 0 ) P+ v / 3}- Since g £ B q , from the definition of the B q norm, 
\g' (L _ ( p+ 1 )/ 2 w) | < || ff || and since $(0) = 0 it follows that \g (L~^ +1 '>/ 2 W ) \ < 
£-(p+ 1 )/ 2 | u; |||g||. Using these bounds, if L > Li we obtain, for all n = 0,1, 2,..., 


\( R L,n9)(w)\ + \(R 0 L,n9)'(w)\ < (l + H + 3(p + 1) VH L (P+1)/2 ||s|| ■ 

-it , 2 

Defining C = sup^ (l + |w| + 3(p + l) _ 1 |w| 2 ) (1 + |'u;| 9 )e 6 <p+ 1 > finishes the proof. I 


The RG analysis involves decomposing the initial condition into two factors: one 
in the direction of the asymptotic fixed point of the RG operator and the other in 
a direction which is irrelevant in the RG sense. That is, when the RG operator is 
applied to the initial data, the irrelevant component is contracted for large L. 


In Lemma 13.41 we decompose f n given by (|30ll and obtain estimates needed later 
to prove Theorem 11.11 We will first assume that, given k £ N, /„ is well defined 
for all n = 0,1 ,..., k, which is guaranteed if ||/ ra || < e n for all n = 0,1 ,..., k — 1 
(cf. Lemma 13.11) . Later, in Lemma 13.61 we prove that if fo is small enough, then 
the sequence {/„} given by (1501) is well defined. Furthermore, notice that from the 
definition of the B q norm and some estimates used in the proof of Lemma 13.31 we 
obtain constants C p , q and K PA such that ||/*|| < C Pyq and, if L > L\, 

\\R%f;\\<K p , q , Vn = 1 , 2 ,... ( 42 ) 

For the next lemmas, we will always refer to the constants C p ^ q and K p q and to C 
and Li given in Lemma 13.31 


Lemma 3.4 (Renormalization Lemma) Given L > L\, suppose f n is defined 
for n = 0, 1, + 1 as specified in (Iddl) . We can then write 

fo = A 0 f;+g 0 , f n+1 = A n+1 R% +1 f;+g n+1 (n = 0,1,..., k) (43) 
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in terms of functions g n G B q , g n ( 0) = 0 and constants A n , K which satisfy 


||ffra+i|| < CL-^ +1 ^ 2 \\g n \\ + iCL n [ p+3_Q ( p+1 )]/ 2 ||/„|| 2 (44) 


and 

\A n+1 - A n \ < C„L n[p+ 3 - Q(p+ 1 )I/ 2 ||/n|| 2 , 


with C n given by $33 1) . 


(45) 


Proof: We first prove (1151) inductively. Define go by fo = A 0 f* + go, with 
Aq = /o(0) and since /*(0) = 1, we have 5o(0) = 0. By hypothesis, /i is well 
defined by Rl,o/o and using representation (1361) and the decomposition above 
for fo we can write /i = AiR^f* + < 71 , where A\ = A 0 + t'o(O) and gi(x) = 
R° L go(x) + L( p+ 1 )/ 2 i/ 0 (-b^ p+1 ^ 2 aO — vo{<S)R° L f*(x). It follows from the definition 
of Rf that tF(R ( lgo)(0) = 0 and tF(R { ff*)( 0) = 1 and, therefore, <?i(0) = 0, which 
proves 0 for n = 0. Now suppose 0 holds for n = 0,..., j — 1, where j < k. We 
will prove that it holds also for n = j. Using (l43l) with n = j— 1, representation (l36l) 
and the semi-group property of the linear RG operator we obtain 

fj+i(x) = AjR 0 L j +1 f*(x) + R° L gj(x ) + ( 46 ) 


Defining 

Aj +1 = Aj + Uj( 0) (47) 

and 

9j+i(x) = R° L gj( x ) + U (p+ 1 )/ 2 ^(i (p+ 1 )/ 2 a:) - ^(0)i?°i+i/*(»), (48) 

we can write dlCH) as fj + 1 = Aj + iR° Lj+1 f* + gj+\. From the induction hypothesis, 
gj( 0) = 0 and therefore, from definition (|48|) . since the Fourier Transforms of Rfgj 
and R^ j+1 f* at the origin are equal, respectively to <jj (0) and /*( 0), we obtain 
5 .7+1 (0) = 0, which proves ( 1451 ) for n = 0, 1 ,..., k. 


Recalling that v n {x) = N n (u)(x,L) and since estimate (1511) holds for all n, using 
definition 0 we obtain 0 for n = 0,1,..., k. After a calculation similar 
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to the one in the proof of Lemma 12.11 we obtain ||Z>+ 1 )/ 2 i/„(Z>+ 1 )/ 2 -)|| < 

L (p+l)q/2 CnL n{p+3-a(p+l)]/2^ n ^2 and using gJJ , 11 ^ (p+1)/2 ^ (p+l)/2.) _ 

^(0)i?°„ +1 /;|| < (L (p+1) « /2 + K p , q )C n L^+ 3 -^P +1 W 2 \\f n \\ 2 . Since g n { 0) = 0 and 
L > Li, from definition (j48|) and Lemma ITTTTTl we obtain 

IlSn+ill < CL~^/ 2 \\g n \\ + (l(p+D?/ 2 + K p ^ C n L n ^+ 3 - a ^ +1 ^ 2 \\f n \\ 2 (49) 

for all n = 0,1,..., k. Now it follows from (14T1) that the constants C n are uniformly 
bounded. In fact, defining 

K = S(L-\) ^2^TI) + l) ( L(p+1)q/2 + K P , q ) S 2 (p 0 ), (50) 

then C n < K for all n and therefore we obtain inequality (l44l) , which ends the proof. 


The estimates obtained in Lemma lXTl are used to prove Theorem 1 1.1 1 in the following 
way: (l45l) guarantees that the sequence (A n ) is convergent and (l44l) is used to 
prove that the component g n gets smaller as we increase n. This is so because 
of our definition of a or, in other words, because the nonlinear perturbation F of 
problem ([3]) is irrelevant. Before we apply Lemma f3.4l we have to prove that the 
initial data of each problem (1?71) is small enough and to do that we will define a 
recursive sequence (G„) such that, for all n, \\f n \\ < G ra ||/o||- In Lemma [3~5l we 
prove that, under a certain condition, this sequence is bounded. Given S £ (0,1), 
let 

L S = max{Li, [2G(1 + Gp,,)] 2 /^} (51) 

and for L > L$, define 

OO 

G = 1 + K p . q L J ' (<5_1)/2 < oo, (52) 

3=0 

Gi = L^ -1 )/ 2 + K Ptq ( 1 + Goll/ll) and G n +i, for n = 1, 2,3,..., by the relation: 

G n+1 = l^-DG+D/ 2 + K p , q ^1 + Goll/H +^C j G 2 L^ +3 ~ a ^y 2 \\ f \\ 
where each Cj , with j = 0,1, 2,..., is given by equation (1351) . with n = j. 
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Lemma 3.5 Let 6 £ (0,1) be such that <5 — 1 > p + 3 — a{p + 1) and let L > Ls, 
where Lg is given by m■ Also, let K and G be given by \50\) and \52 1) . respectively, 
and suppose that f satisfies 

A ' G2 ll/ll < vjhm- ( 53 > 

Then G n+ 1 < G for all n = 0,1,2,... . 


Proof: Since L > 1 and G > 1, it is straightforward from the fact that Co < K and 
from condition (53} that Gi < G. Now, suppose G n+ 1 < G, Vn = 1,2,..., k — 1. 
From the definition of G n +i, using the induction hypothesis and since C n < K,Vn, 


G k+ 1 < i(«-i)(fc+D/2 + Kp q | i + + KG 2 \\f\\ E £p[p+3-a(p+l)]/2 J 

i=1 


Now, from (ESI) and since L > 1 and (5—1 > p + 3— a(p + 1), we obtain 
Gfc+i < 1 + K p ^ q (l + H-1- L < G , which completes the 


proof. 


In Lemma [3.111 we will obtain estimates for the rescaled solutions to IVP (1271) . In 
fact, we will define e > 0 such that, if the initial data / of problem Q is in the 
ball of radius e, then there is a unique global solution to IVP E]). Furthermore, we 
will prove that, under certain hypotheses, the component g n of the initial data /„ 
goes to zero when n —> oo. This fact will be used to obtain the asymptotic behavior 
in Theorem 13.11 Before stating the lemma, we notice that from (l4lT) , if e n is given 
by (135)) and 


a = min 


(2 K)~\ 


2 C q 



3 LP+ 1 \ 

2(p + i)y 



(54) 


then a < e n for all n. In the next Lemma wc will refer to K, Lg, G and er given, 
respectively, by (1501) . (1511) . (1501) and (l5Tll . 


Lemma 3.6 Let L > Lg and S € (0,1) such that 5—1 >p + 3 — a(p + 1). Then, 
there ise> 0 such that, if ||/o| <£, f n given by ISM is well defined for all n> 1, 

||/n||<G n ||/ 0 || (55) 
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and if g n is given by the decomposition 143\) > then, 



(56) 


Proof: The proof is by induction in n. First, define 


e = mm 


jerG” 1 , [2i^G 2 L (1 ”' 5)/2 ]- 1 | 


(57) 


Since G > 1, we have ||/o|| < cr < e 0 and, from Lemmas 13.11 and 13.41 /i is 
well defined by fi = AiR^f* + gi and gi satisfies da with k = 0. Therefore, 
since f 0 = A 0 f* + g 0 , we obtain ||si|| < [G(l + C p , q )L-( p+1 '>/ 2 + if||/o||]||/o||- 
Since L > L$, then, 2G(1 + C PA )L~G +1 )/ 2 < if 15-1 )/ 2 and since G > 1 and 
II/o|| < e, then 2iL||/ 0 || < L (l5_1)/2 . Therefore, ||ffi|| < L (5_1)/2 ||/ 0 ||. Now, using 
decomposition (1431) with k = 0 and the bound (142[) . 


||/i||< (1 + Co\\fo\\)K p , q + l/ 5-1 )/ 2 H/oll = Grll/oll, 


which proves the Theorem for n = 1. Now suppose there exists k > 1 such that (l55l) 
and (1551) hold for all n = 1,2,..., k. We will prove that these estimates hold also 
for n = k + 1. From the induction hypothesis and Lemma 13.51 since \\fo\\ < £, we 
have ||/„|| < G||/ 0 || < e„, Vn = 1, 2,..., k. Therefore, we can apply Lemma [5~H to 
obtain estimate (l44l) with n = k. Then, using (l55l) and (l56l) with n = k, we get: 


C £Mp+3-a(p+l)]/2 


||5fe+i|| <i (5_1)(fc+1)/2 


KGl || /o|| ||/o||. 


L (p+<5)/2 + L (S-l)(k+ 1)/2 


Since G p>g > 0 and L > L^, then CL~^ P+S ^ 2 < 1/2 and since ||/o|| < e, using 
Lemma l3.5l we obtain (1561) with n = k + 1. By Lemma 13.41 fk+i is well dehned and 
can be represented by (l43l) . Therefore, using the triangle inequality and (l42l) . 



Now, since |A 0 | < ||/ 0 || and G„ < K , for all n, applying estimates (l45l) and (l55l) for 
n = 0,1,2,..., k and using Lemma l3.5l we obtain (1551) with n = fc + 1. In particular, 


||/fc+i || < G|| /o|| < e fc+ 1, which ends the proof. 
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We have proved that for |j/ 0 || < e each IVP (1271) has a unique solution u n in Bf n . 
To finish the proof we only need to concatenate these solutions to obtain a unique 
global solution to IVP ((3|). 

We first extend the definition ([5]) of the B^ L> space by considering the space 
S (Ln+1) = {u : R x [L n , L n+1 ] -» R | ||u|| L „+i < 00} 
with the norm ||u|| £ »+i = sup t£[i n £ n+i] ||u(-,i)||. 

Now define {h n } by 

ho = f and h n+1 = L-«(p+ 1 )/ 2 Url ^ £ -n(p+i)/2 a . j L j 

and let Uh n be the solution to IVP (1771) with X n = 0 and initial condition h n . 
Finally, define 

B = {u G B<°°> : ||u - u hn \\ L n + i < ||/i„||,Vn} , (58) 

where (abusing notation) || • || £ n+i denotes the seminorm induced by the obvious 
projection from onto ij( L " +1 ). 

Corollary 3.1 Under the hypotheses of Lemma \3.6\ the IVP has a unique 
solution u € B. In this case, the RG transformation has the “semi-group property” : 
RL*,of = Rl,ti- 1 o • • • o R L1 o R L ,of for all n> 1. 

Proof: From Lemma 1517711 since ||/o|| < e, then ||/„|| < e„ for all n and using 
Lemma 13.11 we obtain the existence and uniqueness of solutions to problems (1271) 
in Bf n . Now define u(x,t) = L~ n( - p+1 ' ) ^ 2 u n (L~ n( - p+1 ^ 2 x, L~ n t) , Vi £ [L n ,L n+1 ] 
and take y = L~ n ( p+1 ' > / 2 x and r = L~ n t. Since u n {y,T) is the unique solution to 
IVP (1771) in Bf n , then u(x,t) is the unique solution to IVP @ in B. To prove the 
semi-group property, it is enough to apply Lemma HO and (1771) . inductively. I 

The previous results are concatenated in Theorem 13.II 
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Theorem 3.1 Under the hypotheses of Lemma \3.6\ there is a constant A such that 


lim 


L n(p+l)/2 u ( L n(p+l)/2^ L nj _ 


= 0 . 


Proof: Since ||/o|| < e, it follows from Corollary [3d] that the IVP J3]) has a unique 
solution u £ B. It follows from the semi-group property and Lemma 13.41 that 
f n = A n R? Ln f* + g n = i"(p+ 1 )/ 2 u (X Tl (p+ 1 )/ 2 x, L n ). Therefore, estimate (l56l) can be 
written as ||/„ - A n R%f *|| < L ti(,5_1)/2 ||/o||. Since C n < I< and ||/ 0 || < e, using 
Lemma 1531 and estimates fl5l) and (1531) . we obtain, for all n = 1, 2, 3,..., 


I^n+l A n | < 


j^n[p+3-a{p+l)\/2 

2LU~ S )/ 2 


ll/oll 


and therefore, ( A n ) is a Cauchy sequence in R. Let A £ 1 be the limit of this 
sequence. Using the triangle inequality, \\L n ^ p+1 ^ 2 u(L n ^ p+1 ^ 2 x, L n ) — Af*\\ < 
\\L^+V/ 2 u(L n(p+1)/2 X ,L n ) - A n R%f ;|| + \A\\\R° Ln f* - f ;|| + |T„ - A\\\R%f;\\, 
which, from Lemma 13.21 and (1421) can be upper bounded by 


\A\M 


r(L n ) 

£ j n[p+3-a(p+l)]/2 

L n (P+ !) 

' 2L ( ' S 1 ~ S)/2 (l - /4 p+3_a(p+1)I/2 ) 


K P , q \\fol (59) 


T«(l-<5)/2 

Then it is enough to take the limit when n —> oo. I 

Theorem o now follows from estimate (1591) as we explain below. 

Proof of Theorem ll.lt We have proved that (j4j) is valid when the initial data 
/ = /o is sufficiently small and t = L n (n = 1,2,...), for L > Lg. In fact, it 
follows from (1551) that, if t = L n , then || Vt p+1 u(Vt p+1 x : t) — Af*\\ < ^U -1 )/ 2 1|/ 0 || + 
|T|M|f-(P +1 )r(t)| + fb+a-aCP+ill/s^U- 1 )/ 2 2 (l - 4 P+3_ “ (p+1)]/2 ) 1 K p ^\\f 0 \\. 

We can extend this bound to t = rL n , with r £ [1 ,L\ and L > Lg by replacing 
everywhere L by r l / n L. Therefore, since the constants in (159)) do not depend on 
the particular value of L > Lg considered, the inequality above holds for all t > Lg. 
Taking the limit t —* oo finishes the proof. I 


To conclude, we remark that we do not have an explicit expression for the limit A 
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of the sequence A n . However, it should be clear that 


A = lim / u(x, t) dx. 
t^oo J 

Also, we point out that the RG approach has been used to obtain higher order 
corrections to the asymptotic behavior, see 00)- 
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